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Abstract 

The existence of a pseudogap above the critical temperature has been widely used 
to explain the anomalous behaviour of the normal state of high-temperature super- 
conductors. In two dimensions the existence of a pseudogap phase has already been 
demonstrated in a simple model. It can now be shown that the pseudogap phase 
persists even for the more realistic case where coherent interlayer tunneling is taken 
into account. The effective anisotropy is surprisingly large and even increases with 
increasing carrier density. 
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1 Introduction 



The anomalous behaviour of the normal state of high-temperature superconductors (HTSC) 
|1|, 0] (including the behaviour of the spin susceptibility, resistivity, specific heat and photo- 
emission spectra) has been recently interpreted in terms of the formation of a pseudogap 
above the critical temperature, Tc @, 

The formation of a pseudogap phase above has been explicitly demonstrated in a 
model non-relativistic 2D fermi-system The work is based on the peculiarities of 

the Berezinskii-Kosterlitz-Thouless (BKT) phase formation (see also 0), which is a two 
stage process. For a 2D system one must rewrite the order parameter $(x), where x = 
r, r denotes the position and imaginary time, in terms of its modulus p{x) and its phase 
6{x) i.e. $(a;) = p{x) exp[—i9{x)]. This was originally stated by Witten in the context of 
2D quantum field theory [Q. It is clearly impossible to obtain $ = {^{x)) 7^ at finite 
T since this corresponds to the formation of homogeneous long-range (superconducting) 
order which is forbidden by the Coleman-Mermin-Wagner-Hohenberg theorem. However it 
is possible to obtain p = {p{x)) 7^ but at the same time $ = p{exp[—i9{x)]) = due 
to random fluctuations in the phase 6{x). We stress that p 7^ does not imply long-range 
superconducting order (which is destroyed by the phase fluctuations) and is therefore not in 
contradiction with the above-mentioned theorem. 

Thus one has three regions in the 2D phase diagram [^]. The first is the superconducting 
(BKT) phase, where p 7^ 0. In this region there is algebraic order and a power law decay 
of the correlations. The second is the pseudogap phase where p is still non-zero but the 
correlations decay exponentially. The third is the normal (Fermi-liquid) phase where p = 0. 
Note that $ = everywhere. The unusual properties of the second region, which lies be- 
tween the superconducting and normal phases, have previously been used in P| to explain 
the pseudogap behaviour in HTSC. For example in the mean-field calculation of the para- 
magnetic susceptibility P] the parameter p plays the role of the energy gap A in the theory 
of ordinary superconductors. Thus in this calculation one has the opening of an energy gap 
p (or equivalently a lowering of the density of states) above the superconducting transition 
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temperature. 

The description of the phase fluctuations and the BKT transition closely resembles that 
given by Emery and Kivelson However in their phenomenological approach the field p{x) 
does not appear while in the present microscopic approach based on ^ it appears rather 
naturally. 

Although the description in terms of modulus and phase variables is essential from a 
mathematical point of view, the physical implications for the theory of superconductivity 
have not been fully understood. Previous work |^ was at zero temperature but the applica- 
tion to condensed matter requires the extension to finite temperature where p is a function 
of the temperature T. It makes sense therefore to define the temperature Tp at which p be- 
comes zero. This temperature is then interpreted as the temperature at which the pseudogap 
opens. Since in p, 0] and in the present work p{x) has only been treated in the mean-field 
approximation i.e. one has neglected the fluctuations in both p{x) and 6{x), a second-order 
phase transition was obtained. However it is well known experimentally that the formation 
of the pseudogap phase does not display any sharp transition. It can be argued however that 
the fluctuations may convert the obtained sharp transition to a crossover 

An important question is whether the pseudogap phase (PP) forms in real HTSC, which 
are only quasi-2D systems. Another related question is whether the pseudogap phase, if 
present, remains large enough to explain the experimentally observed anomalies even when 
interlayer tunneling is included. 

This paper will study PP formation for a quasi-2D system with arbitrary carrier density. 
Note that this problem is further complicated by the formation in quasi-2D systems of the 
ordinary bulk superconducting phase with homogeneous long-range order (LRO) at critical 
temperature Tc (see for example [|l3])- Since the BKT phase formation corresponds to 
the formation of two-dimensional order while true LRO is three dimensional one expects 
that for weakly coupled layers Tc < T^kt^ where Tgj^rp denotes the critical temperature for 
the BKT transition in the quasi-2D system. Tgj^rp represents the temperature at which the 
system becomes superconducting in the layers (planes) while is the critical temperature for 
bulk superconductivity. Thus Tgj^-p is the maximum temperature at which superconducting 
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behaviour is present. The temperature range Tc < T < T^^^t corresponds to a crossover 
region where the effect of the interlayer tunnehng is insufficient to produce three dimensional 
behaviour and the behaviour remains BKT-hke. There is in fact some experimental evidence 
for the existence of this crossover region (~ IK) in optimally doped YBa2Cu307_5 fTT ]. 



The size of this region is obviously critically dependent on the anisotropy, the carrier 
density and the precise form of the interlayer tunneling. At very low densities (i.e. in 
the Bose limit of isolated pairs) and for physically reasonable values for the parameters, it 
has been shown for quasi-2D system |12|, ^ that long-range order is only established well 
below Tbkt- However these densities are not realised physically. It can also be shown that 



approaches the BCS critical temperature asymptotically in the high-density limit so that 
both this region and the pseudogap phase must vanish asymptotically. At the intermediate 
carrier densities found in HTSC we find that the pseudogap phase remains large enough 
to explain the observed anomalies. However the experimentally observed superconducting 



transition is practically three dimensional |TT| so that one expects Tc<T^^r^ at these densities 
and this is the subject of current investigations. 

For these reasons we only calculate the temperatures Tgj^rp and Tp as functions of the 
carrier density Uf to establish the boundaries of the PP. Here T^j^rp denotes the critical 
temperature for the BKT transition in the quasi-2D system. 

We show that the value of Tp is practically identical in the 2D and the quasi-2D systems, 
while for realistic model parameters Tgj^rp > Tbkt- More importantly, even at relatively 
low anisotropy (for example that observed in the compound YBa2Cu307_5), the difference 
between Tg^^r^ and Tbkt is too small to destroy the PP at all carrier densities. In addition the 
difference at high doping levels can be shown to tend to zero logarithmically. This enables 
one to seriously consider the pseudogap phase explanation claimed in reference |^ even in 
the quasi-2D case. 
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2 Model and Formalism 



The nature of the interplane tunnehng in HTSC is not yet well established and several 
different models exist. Here we choose the simplest possible Hamiltonian density, often 



employed for studying HTSC [|T^, |1G 
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H 15- cos{idVz) + IJ' 



ip^{x) - Vi)\{x)i)\{x)i)i{x)i)^{x), 



2m± nizd? 

where x = T,rj_,rz (with r_|_ being a 2D vector); is a fermion field, a =t, i is the 

spin variable; m± is the effective carrier mass in the planes (for example CUO2 planes); 
is an effective mass in the 2;-direction; d is the interlayer distance; V is an effective local 
attraction constant; /i is a chemical potential which fixes the carrier density nj; and we take 
h = kB = I. 

The Hamiltonian proposed proves to be very convenient for the study of fluctuation 
stabilization by weak 3D one-particle inter-plane tunneling. We have omitted in (|I|) the 
two particle (Josephson) tunneling considering it to be less important than the one-particle 
coherent tunneling already included. There can be situations where Josephson tunneling is 
more important. In fact some authors consider the most important mechanism for HTSC to 
be the incoherent inter-plane hopping (through, for instance, the impurity (localized) states 
or due to the assistance of phonons). We will not however consider Josephson tunneling 
here. We do however take into account the layered structure of HTSC which is a vitally 
important extension to the 2D models usually considered. 

It is significant that the large anisotropy in the conductivity cannot be identified with 
the corresponding anisotropy in the effective masses rriz and m±. In particular, HTSC with 
rather large anisotropy in the z-direction do not display the usual metal behaviour at low 



temperatures [17|. However this semiconducting behaviour is not directly related to the 



pseudogap phenomena and the Hamiltonian (|l]) may be used to study the qualitative 



features of pseudogap opening. 

The Hubbard- Stratonovich method was applied to study the system described by (||). In 
this method the statistical sum Z{v,fi,T) is given as a functional integral over the Fermi- 
fields (Nambu spinors) and the auxiliary field $(x) = In contrast to the usual method 



where one calculates Z in terms of the and variables, the parameterisation = 

p{x) exp \—id{x)\ should be used ^ (see also [jl8|, |19[). In addition to this reparameterisation 
one must make the replacement V^o-l^;) = Xaip^^^^V This representation splits the 

charged fermion field ^/'o-(a;) into a neutral fermion field Xa{p^) and a charged boson field part 
exp [i^(x)/2]. This resembles the spinous and holons in Anderson's approach. 

This particular choice of parameterisation ensures that is single- valued with period 
2tx. As a result one obtains 



Z{v,p,T) = J pVpVOexp [-pn{v, p,T, p{x),d9{x))], 



(2) 



where 



Q = - f dr f drp^ - TTr LnG"^ (3) 
V Jo J 

is the one-loop effective action which now depends on the modulus-phase variables. The 
action (^) is expressed in terms of the Green function G of the initial (charged) fermions 
which now has the following operator form G^^ = Q^^ — S where 



G-M^)] = -idr + r^ 

i 



' H cos{id\/z) + p 



2m 



m 



Tip{x); 



(4) 



'idre , {v^ey , (v,g)^ 

1 1 cosuctV, 



iv^.e ivle 



iVzO sm{id'Vz 
2m,d 



(5) 



4m_L Aniz 2m± 
Here Q is the Green function for the neutral fermions. 

Note that in S we have omitted higher order terms in VzO but in order to keep all 
relevant terms in the expansion of sm{idVz) the necessary resummation was done. Since 
the low-energy dynamics in the phases in which p 7^ is determined by the long-wavelength 
fiuctuations of 0{x), only the lowest order derivatives of the phase need be retained in what 
follows. This gives the one-loop effective action as 



n nUv, T, p, 09) + n^^{v, p, t, p) 



where 



^kiniv, Ai, T, p, 86) = TTr ^ -(^S) 



n=l 



n 



(6) 



(7) 



p — const 



and 



'pot 



[v, /i, T, p)= j dvp^ - TTr Ln^ ^ 



Given the representation (1) one can obtain the full set of equations for Tgj^rp, p(Tg 



p — const 

IS for T^KT, P(^BKTJ 

and /i(7Bj^rp) at given ep. While the equation for Tgj^rp only depends on the kinetic part ([^) 
of the effective action, the equations for p(Tgj^rp) and y^iT^^:^ can to a good approximation 
be obtained using the mean field potential @. In the phase where p 7^ the mean-field 
approximation for p describes the system well due to the nonperturbative character of the 
Hubbard-Stratonovich method. 

We note that the expression for the potential (§) in terms of p^ is identical to the mean- 
field potential in the BCS approximation but with |$p replaced by p^ [Q. Thus Tp, the 
temperature at which p = 0, is in this approximation identical to the BCS mean field 
temperature T*^^. However, although Tp = T^'^^ in the mean-field approximation for p(x), 
the two temperatures have a very different basis, both mathematically and physically. This 
becomes evident when one includes the fluctuations. Not only does Tp remain finite in 
two dimensions due to the structure of the perturbation theory in the new modulus-phase 
variables, but it is also bounded below by Tbkt- On the other hand T^^p approximates the 
temperature where $ becomes non-zero (onset of long-range order) which is zero in two 
dimensions. 

The modulus-phase representation introduced here is a good tool to consider different 
types of short range order. The success of the BCS approximation is related to the fact that 
in 3D system with large carrier density short and long range order set in simultaneously 
because the fluctuations do not change the situation drastically. 



3 The Berezinskii-Kosterlitz-Thouless transition in quasi- 
2D theory 

If the model under consideration reduced to some known model describing the BKT phase 
transition, we could easily write the equation for Tbkt- Indeed, in the lowest orders the 
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kinetic term coincides with classical spin quasi-2D XY-model (see also |Ty]) which 
has the following continuum Hamiltonian 

H = ^J dr{J[V^e{r)r + UVAr)?}. (9) 

Here J and are constant coefficients (J^ <^ J). Unfortunately the temperature for the 
BKT transition in the quasi-2D case is not as well investigated as in the pure 2D case. In 
fact only in the highly anisotropic case, a = Jz/ J <^ 1, when the vortex ring excitations are 
irrelevant has the transition temperature been derived [^. In this limit the temperature 
^BKT f*^^ BKT transition in the quasi-2D system is close to that in the pure 2D case and 
determined by the equation 

TT r Rtt 

(10) 



BKT — ~ 



2 



»7r 
m a 



This equation was given in |^ and employed to calculate Tgj^-p for the relativistic quasi-2D 



four- Fermi theory |T0[. The equation was derived using the renormalization group technique, 
which takes into account the non-single-valuedness of the phase 6. Thus, the fluctuations 
of the phase are taken into account at a higher approximation than Gaussian. Below the 
temperature T^x^t ^^e correlation function 

while above T^^^r ^^is correlator decreases exponentially. This is the BKT transition in the 
classical quasi-2D XY model. The small correction to the unit in the brackets of Eq. (|TUp 
corresponds to the influence of the third direction. We note that the equation is only correct 
in the limit a <^ 1. From a physical point of view it is evident that the BKT phase cannot 
form for a ~ 1, which corresponds to the 3D limit. Moreover, if the temperature T^. of the 
LRO formation approaches Tgj^rp from below, the BKT phase will not form and one will have 
a superconducting transition directly into a phase with LRO. This may well correspond to 
the experimental situation and we hope to study this question in detail in our future work. 

To expand fi^m up to ~ (V^)^, it is sufficient to consider only the terms with n = 1, 2 in 
the expansion (|^). The method is the same as that in I, P, and gives 
T r'^ 



2d Jo 



I dr J dr [j{fi,T,p{fi,T)){V±e)' + Ufi,T,p{fi,T)){Vzey+ 



K{^^, T, p(/i, T)){dreY + n^(/i, T, p(/i, T))tdre 



(12) 



where 



J {fx, T, p) 



d 



Am 



-nFilx,T, p) 



dt 



,„x + ^/2r + «,cost/(2r) ^j3j 



A* W cost / 2 



4T^ 



Jz{p^,T, p) 



m± 1 



Am, (27r)2 



J dt cos t |w cos t + \l {w cos t + /i)^ + p2 



+ 2rin 



1 + exp 



^ti'cost + pY + 



w I dt sin^ t 
Jo 



J£ w cost 

2T 2'i' 



cosh ^ / x'^ 



AT' 



K{p,T, p) 



m_|_ 



W cos t + 



(14) 



wcost + jiY + , 
tanh^^^ r.^-^ ^1 (15) 



(47r)2 7o ^ ' y^(^i;cost + p)2 + p2 " 2T 
and where we have introduced the bandwidth in the ^-direction w = (m-,d'^)~^. Here 



nF{p,T, p) 



m± 
(27r)2d 



+ 2Tln 



J (it |/i + {w cost + /i)^ + p2 



1 + exp 



w cost + p)2 + p2 



(16) 



takes the form of a fermi-quasiparticle density. Note that J(p, T, p = 0) = J2(p, T, p = 0) = 
0. This property of the phase stiffness is present in the 2D model where it imphes that p 7^ 
at Tbkt and thus Tp > Tbkt- The continued presence of this property in the quasi-2D case 
imphes that the essentially nontrivial cosine dispersion law for motion in the third direction 
is correctly treated. 

Now we discuss the features of the BKT transition in our model. In contrast to the 
Hamiltonian (y) the field 6 is also depend on the imaginary time r, and therefore by Fourier 
decomposition, one can write 



6{t, r) = ^ exp(i27rnTr)0„(r.) 



:i7) 
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From Eq. (|T2|) we can see that the nonzero mode 6n{r) n ^ has a mass ~ {2iTnT)'^K, 
and only the massless component ^o(r) is relevant in the low-energy region. In terms of the 
9n fields, the effective classical action is expressed by 



Recalling that 6 is an angular variable, we see that the ^o-P^rt of the effective classical action 



(p^) is nothing but the Hamiltonian of the quasi-2D XY model (^. This makes it possible 
to write the equation for Tgj^rp: 

Stt 



71" 

^BKT = ^BKT) P) 



1 + 



(19) 



where a{fi,T,p) = J;,{fi,T, p) / J{fi,T, p) is a function of p, T and p{p,T). Recall that 
equation (|l^) is only correct in the limit a <^ 1. Thus after the calculation of Tgj^rp one must 
check that the condition a -C 1 is satisfied. 

6o has non-trivial dynamics described by the correlation function 

^gieo{r)g-ieo(o)^ ^ ^T/4.j ^ ^ ^20) 

which is identical to ([TT|) but with 6 replaced by 6q. However the gauge- invariant correlation 
function, ($*(r = 0, r)$(0)), has an additional factor from the massive modes n 7^ and, 
regardless of the correlator for 6^0 (|20|) , one has 



($*(r = 0, r)<l>(0)) -^0 r ^ 00 (21) 

for any T > (see the details in [^]) which is consistent with the assumption that LRO is 
absent. 

Although mathematically the problem reduces to a known problem, the analogy is in- 
complete. Indeed, in the XY model the vector (spin) subject to ordering is assumed to be a 
unit vector with no dependence on T. In our model this is not the case, and a self-consistent 
calculation of Tbkt as a function of nj requires additional equations for p and p, which 
together with ([T9|) form a complete system. 
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Using the definition (H), one obtains (see e.g. 0) 

^2 " dk 



apot{v,fi,T, p) 



V 



2T In cosh ' \L - ^(k) 



(22) 



(27r)3 1 2T 

where .^(k) = £:(k) — fi with £:(k) = k^/2m_L — w cos kzd. Then the missing equations are the 
condition dQpot{p)/dp = and the equahty v~^dQpot/dp = —rif, which fixes n/ : 

dk 



1 
V 



1 



tanh 



e(k)+p2 

2T 



(23) 



nF{p,T,p) =nf. (24) 

The equations ( P^D and ( P^ obtained above comprise a self-consistent system for determining 
the modulus p of the order parameter and the chemical potential p in the mean-field approx- 
imation for fixed T and Uf. In the phase where p ^ the mean-field approximation for p 
describes the system well due to the nonperturbative character of the Hubbard-Stratonovich 
method and the character of the perturbation theory in the modulus-phase variables. 

To simplify the problem one can take the limit — oo in the expressions (|T3|), (|16D 
assuming that mz/m± ^ 1 and w = m~^d^'^ <^ ^bkt- This is indeed the case in real 
HTSC compounds e.g. for ~ lO^mg and d = lOA the value of / {iiizd^kB) ~ lOK is far 
less than the critical temperature. This simplification is useful but not essential since the 
numerical results are practically unchanged by this approximation. 

The energy of two-particle bound states in vacuum = — 2W^exp(— 47r(i/m_L\^) (see e.g. 
||22|| ) where W is the bandwidth in the plane, is more convenient than the four-fermi constant 
V . For example, one can easily take the limits W oo and — in the equation (p3|), 
which after this renormalization becomes (in the limit oo) 

\eb\ _ /-^ . 1 



In 



V/i^ + p2 - p 



2 / du- 



'^2 + 



exp \l + 



(25) 



Thus, in practice, we will solve numerically the system of equations (p!9D, (p5|) and ( p4D to 
study T^j^rp as a function of n/. 

Setting p = in the equations (p3D and (p^ , we arrive (in the same approximation) at 
the equations for the critical temperature Tp above which p = and the corresponding value 
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of /i: 



ln^- = -/ rfM 7 = 1-781, 26 

Ip n Jo u 



1 + exp ( 



T , 



eF, (27) 



where ei? = Tmfd/m± is the Fermi energy H. Note once more that these equations coincide 
with the system which determines the mean-field temperature J'J2D)mf /i(Tj^^^*^^) (see 
\ T\ ^). This is evidently related to the mean- field approximation used here and the limit 
oo. 

Certainly in the simplest Landau theory one appears to have a phase transition since 
p takes on a non-zero value only below Tp. In fact the parameter p describes a gap only 
in the spectrum of the neutral fermion field Xcr{x)- In describing the charged (physical) 
fermion field il)fj{x)., p only appears in conjunction with 9 in every correlation function. At 
zero temperature in closely related four-Fermi models P, ^ the pole in the neutral fermion 
Green's function associated with the gap in the neutral fermion spectrum is converted by 
the phase fiuctuations into a branch cut in the Green's function for charged fermions and we 
strongly suspect that similar behaviour is present at finite temperature. 

In the present approximation the neutral fermion Green's function ^(k, u;) (^) is identical 
to the BCS Green's function but with energy gap p rather than A. As such the spectral 
density of the neutral fermion Green's function, is the sum of two delta functions centered 



on the isolated poles of ^, = ±-E(k) where E{k) = ^^^(k) + and one has zero density 
of states inside the gap p |2^. In the case of a branch cut (as postulated for the charged 
fermion Green's function) the spectral density is non-zero for uj in the entire range — £'(k) < 
u < -E(k). Thus the spectral density is smeared and one expects a corresponding smearing 
of the gap. An illustrative example where a branch cut does lead to pseudogap behaviour is 
given in 
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We believe that (see also discussion in 0) the expected smearing of the neutral fermion 
gap in the present model may well describe the observed pseudogap behaviour. Furthermore 
the fiuctuations in 6 are expected to convert the sharp phase transition to a smooth crossover. 
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4 Results and Conclusions 



The analytical and numerical investigation of the systems (|T^), (P^D, ( [2^ ) and (PB|), ( P7[ ) 
yield the following results. 

One can show that at large densities {ep \£b\ T) when /i ^ 



J{eF,T,p) = -ep 



1 

1 / dx 

2 J—oo 



cosh^ Wx^ + — 



4T^ 



(2^ 



and 



1 

27r 



-to 



dx- 



cosh^ W + 



4T' 



(29) 



which yields 



=-r = — 



(30) 



which is independent of both p and T. Thus even for a modest anisotropy in the initial 
fermion Hamiltonian one expects to obtain a large anisotropy (or equivalently a low value 
of a) in the effective quasi-2D XY model Hamiltonian. Indeed at all densities we find small 
values for a. Surprisingly a even decreases as the carrier density increases. For this reason 
one can indeed use the equation (|^) to determine Tgj^rp. Therefore one expects the results 
obtained for the pure 2D model to persist to the quasi-2D case. 

Indeed for modest carrier densities (ep'<10^|e;,| or the underdoped case) the PP (see 
Fig. 1) in the present model is roughly equal in size to the superconducting region i.e. (Tp — 
^bkt)/^bkt ~ 1- This allows us to believe that in spite of the oversimplified character of the 
model proposed in ||^ this approach may explain some of the observed normal state anomalies 
of HTSC For these densities we argue that the temperature for true LRO is probably 
roughly equal to Tgj^rp. For very low densities we find a large region between Tc and Tgj^rp 
where the system is superconducting but one has two-dimensional and not three-dimensional 
order. It is difficult to say whether such behaviour can be observed experimentally because 
the densities are so low that the Fermi surface is absent [IBl 1131 . 
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The model considered here is also a simple one. Certainly at sufficiently large densities 
the BKT and PP region will disappear due to a direct transition to the state with long-range 
order particularly for low anisotropy. For example, in the case of an indirect interaction in 



2D, it has already been shown [13[ that the PP region only exists at low carrier density. 

The other results obtained in |Q remain valid. In particular one obtains a linear depen- 
dence of the critical superconducting temperature (T^kt) ^^e carrier density over a wide 
range of densities as is seen in experiment. 

For instance the kink in /i at T = Tp, (discussed in reference [^) occurs at the NP-PP 
boundary or before superconductivity appears. The ratio 2A/Tgj^rp is also always greater 
than the standard BCS value The concentration behaviour of this ratio is consistent 
with experiment, namely it decreases with increasing ep- 

It is also interesting to note that the qualitative difference between the temperature 
dependence of the spin and charge correlations in the normal state of the 2D attractive 
Hubbard model which was found in |2^ acquires a natural explanation in the framework of 
the approach used here. Indeed, while the neutral (with spin, but chargeless) fermions have 
the gap p in the PP, this gap should be smeared out for the initial charged fermions due to 
the phase fluctuations. 

To summarize, the presence of interlayer tunneling does not destroy the pseudogap phase 
above the critical superconducting temperature since the effective anisotropy (see (0)) is 
far larger than the simple estimate of m^/m^. 
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Figure 1: T^kt^ ^bkt (dots) and Tp versus the noninteracting fermion density. The regions 
of the normal phase (NP), pscudogap phase (PP) and BKT phase are indicated. Note that 
in the approximation used the value of Tp is the same for the 2D and quasi-2D models. We 
assumed that mz/m± — 100 and {mz(P\eb\)~^ — 0.1. 
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